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ABSTRACT 

Using the gamma-ray bursts simultaneously detected by Swift /BAT and 
Fermi / GBM we performed a joint spectral and temporal analysis of the prompt emis- 
sion data and confirm the rough correlation between the BAT-band photon index 
pBAT g^j^j-j ^j^g peak spectral energy Ep^ak- With the redshift known sub-sample, we 
derived the isotropic gamma-ray energy -E-y^iso and also confirm the -E-y^so — -£'pcak,rost 
relation, with a larger scatter than the Amati sample but consistent with GBM team 
analyses. We also compare the Tgg values derived in the GBM band with those de- 
rived in the BAT band and find that for long GRBs the BAT Tgo is usually longer 
than the GBM Tgo, while for short GRBs the trend reverses. This is consistent with 
the soft/hard nature of long/short GRBs and suggests the importance of an energy- 
dependent temporal analysis of GRBs. 
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1 INTRODUCTION 

Since the launch of the Fermi satellite in 2008, a significant 
number of gamma-ray bursts (GRBs) have been observed by 
both the Fermi Ga mma-ray Burst Mon itor (GBM, 8 keV- 
40 MeV) and Swift (iGehrels et al. 20041 ). Since the energy 
band of Swift BAT is narrow (15-150 keV) and the peak 
of th e spectrum is often o utside the observed energy win- 
dow (|Sakamoto et al. 20091 ). the joint BAT-GBM sample is 
valuable for Swift science by providing important informa- 
tion about the prompt emission, especially the character- 
istic peak energy of the i/Fi, spectrum, -Epcak- Using the 
Fermi /GBM- Swift /HAT joint sample we performed a two- 
part analysis on the prompt emission data. First, a joint 
spectral analysis where we fit the data from both detec- 
tors and investigated several empirical correlations. Second, 
a temporal analysis where we rigorously derived Tgo in the 
8-1000 keV band and compared it to the published BAT 
values. 

A variety of empirical relati ons have been 
discussed i n the literature (lAmati et al. 2002 ; 
Amati 20031: ISakamoto et al. 20041 : ISakamoto et al. 200e : 
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Ghirlanda et al. 2004; 
Yonotoku ot al. 2004), many of which utilize iJpoak- 
Since -Epeak cannot be well measured for most Swift GRBs, 
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efforts have previously been made to look for indicators of 
Epcak based on the available observed quantities. In partic- 
ular, the effective photon power law index in the BAT band, 
pBAT^ has been found to be broadly correlated with -Epeak- 
Two versions of this relation are found in the literature: one 
presented in Sakamoto et al. (2009) using simulations of 
BAT spectra and Ep^ak measurements from broad-band de- 
tectors (e.g. Konus-Wind), and another presented in Zhang 
et al. (2007b) with Ep^ak estimates based on the hardness 
ratio in the BAT band itself (Zhang et al. 2007a; see also 
Cui et al. 2005). In §3.2 we re-calibrate this correlation 
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with the GBM-BAT joint sample. The -E.^.iso 

relation (or Amati relation; lAmati et al. 2003 : lAmati 20031 ) 
is one of the most widely studied, and hotly debated, em- 
pirical correlations that may connect to GRB physics (e 
iBand fc Preece 20051 : iNakar fc Piran 20051 : Ikocevski 201ll 
In §3.3 we test the Amati relation with the redshift-known 
sub-sample of the GBM-BAT joint sample. 

Finally, traditionally long and short GRBs were defined 
in the BATSE band with a separation of about 2 seconds 
in the observed frame (Kouveliotou et al. 1993). When ap- 
plied to Swift GRBs and combined with multi-wavelength 
afterglow data, this definition creates confusion when used 
for GRB classification (e.g. Zhang et al. 2007b, 2009). Since 
the energy band of GBM is similar to that of BATSE, a 
comparison between the measured Tgo durations in the two 
detectors is of great interest and is performed in §4. 
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2 SAMPLE AND DATA REDUCTION 

We conducted our analysis on the sample of GRBs that 
were observed simultaneously with the Swift /BAT and 
_Femij/GBM detectors. The sample was chosen from bursts 
that were consistent both temporally (similar trigger times 
in Swtft MET) and in sky placement (RA/DEC) between 
June 2008 and May 2011. Our final sample consists of 75 
bursts whose spectra were fit individually to check consis- 
tency with the literature and then fit jointly. Some bursts 
were removed because they were a Swift ground detec- 
tion (e.g. 100427A; partial data) or a false detection (e.g. 
080822B which has been classified as an SGR flare). 

The Fermi-GBM GRB data was downloaded from 
the online Fermi GRB burst catalog found on the NASA 
HEASARC websittQ and the time-integrated spectra were 
fit with the GBM software RMfitQ. The Swift-B AT data 
was downloaded from the Swift data table, also lo- 
cated on the NASA HEASARC websitfl The resulting 
data was processed using the typical Swift-BAT products 
(batgrbproduct) after updating to the current calibration 
database (bateconvert) to create the burst lightcurves and 
spectra. 



3 SPECTRAL ANALYSIS 
3.1 Results 

We performed spectral fits for the GBM GRBs with RMf it 
using the built-in power law (PL), power law with exponen- 
tial cutoff (GPL, 100 ke V normalization energy) and Band 
function (|Band et al. 1993 ) models over the entire burst 
temporal and spectral (8 keV-40 MeV) interval in order to 
determine the best fit model. The goodness of the fits are 
measured with the G-statistic. If the G-stat was comparable 
for two models, the simpler model (i.e. the one with less pa- 
rameters) was chosen as the better fit to compensate for the 
added degree of freedom. On many occasions, the decrease 
in the statistic was greater than 30-50 for a PL versus a GPL 
fit, for example, but not significantly reduced when compar- 
ing the GPL and Band function fit. Evidence of this can be 
seen in the majority of GPL fits in the GBM data summa- 
rized in Tables 1 and 2. The spectral regime chosen for the 
NaT detectors was roughly 8-900 keV and 200 keV-40 MeV 
for the BGO detectors. Our spectral fit results are generally 
consistent with the literature (e.g. Bissaldi et al. 2011; Nava 
et al. 2011; individual burst GGN reports). 

For BAT spectral fitting, it is important to make sure 
that the the satellite was not slewing during data acqui- 
sition. If Swift was slewing while capturing a spectrum, a 
new weighte d response file was cr eated to compensate for 
the motion (|Sakamoto et al. 201 ll ). The response matrices 
are created in 5 second intervals and weighted to created 
a final response file that is used with the total spectrum. 
Each time-integrated spectrum was fit over the full BAT en- 
ergy range of 15-150 keV using the same three models listed 
above. If a model showed improvement greater than Ax^ > 
6 per added degree of freedom then that model was deemed 

^ http:/ /heasarc. gsfc.nasa.gov/W3Browsc/fcrmi/fermigbrst. html 

^ fermi .gsfc. nasa.gov / ssc /data / analysis / 

^ http://hcasarc.gsfc.nasa.gov/cgi-bin/W3Browsc/swift.pl 



a better fit. Our spectral fitting for the BAT data and joint 
spectral fits for the GBM and BAT observed spectra are per- 
formed with Xspec. We re-bin the spectra with the criterion 
that the number of photon counts per bin is greater than 
20 by using the HEAsoft tool grppha if the burst was suf- 
ficiently short or weak that this procedure was warranted. 
In order to address the calibration between the BAT and 
the GBM Nal/BGO detectors, each spectrum was fit with a 
varying calibration constant. The BAT calibration was fixed 
to 1.0 and the Nal and BGO constants were left as free pa- 
rameters to be fit. About half of the bursts did not produce 
reliable fits to the BGO detector calibration constant, likely 
due to the sometimes low signal in the BGO detector. For 
these bursts the BGO constant was fixed to the mean value 
of the distribution of BGO calibration constants (1.3) and 
the spectra were refit. The a-nd spectral parameter errors 
were generally reduced when compared to the fits where all 
the calibration constants were forced to be equal. 

The best-fit spectral models of all the 75 bursts in the 
joint GBM-BAT sample are presented in Table 1 (simple 
powerlaw model) and Table 2 (cutoff powerlaw or Band 
models). The distributions of iJpeak and a of the sample are 
shown in Fig. [T] The results are generally consistent with 
previous results (e.g. IPreece et al. 200d : I Zhang et al. 201 ll : 
iNava et al. 20111 ). The Speak distribution is roughly log- 
normal but is somewhat wider than that of the bright 
BATSE GRB sample (Preece et al. 2000). Our sample in- 
cludes GRBs in a broader Eiso range, broadening the ob- 
served -Bpcak range. The low-energy spectral index a peaks 
around 1 as suggested by previous works. We do not present 
the distribution of the high energy photon index /3 due to 
the small sample size (most spectra can be fit by a cut- 
off powerlaw model). Although they form a limited part of 
our sample the distribution for short bursts show behavior 
consistent with previous samples, namely a harder median 
Epoak and similar or slightly shallower pre-break slope. 

Since the GBM energy band covers that of BAT, we also 
derive iJpcak with the GBM data only and compare the val- 
ues with the joint fit values. Figure [2] shows the relationship 
between the separately derived GBM and joint sample Speak 
values with a line of slope 1 drawn for reference, showing the 
generally small differences between the samples. 

3.2 r'^'^'^ - Ep relation 

With our spectral analysis results, we examine the empirical 
relations between the PL slope of the BAT spectrum, F^*^, 
and -Epeak proposed by Sakamoto et al. (2009) and Zhang et 
al. (2007b): 

log Speak = 3.258 - 0.829r''^'^ (1) 
and 

log Speak = (2.76 ± 0.07) - (3.61 ± 0.26) logr^*^,(2) 

respectiveljjf]. Both Zhang et al. (2007a, b) and Sakamoto et 
al. (2009) indicate that the relation is valid between slopes 

^ Sakamoto ct al. (2009) derive Icr error levels which are based on 
the value of T (see their Table 1): Lower Umit: -20.684 -|- 43.646r 
- 26.89ir2 + S.lSSr^ Upper limit: -5.198 + 16.568r - lO.eSOF^ 
+ 2.034r3 
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of roughly 1.2 < r^^^ < 2.3. Here we aim to present the 
entire data set of the joint sample and see how the relation 
fares over the entire observable sample. It is also interesting 
to explore whether using the joint spectra is consistent with 
the derivation of -Epeak from using only the GBM sample 
and a motivation for investigating the combined constraints 
of the two detectors. 

Using only the 53 joint sample bursts that had BAT 
spectra best fit with a PL function and the joint spectra 
Epcak, regardless of whether it is from a CPL or Band func- 
tion, we derive relations that are steeper than the values 
above (Figure [3]). The ordinary least squares method gives 
shallower slopes more consistent with the values of Sakamoto 
(2009) and the bisector fits are steeper and more consistent 
with the values presented in Zhang et al. (2007b). Below we 
summarize the bisector fits, assuming these produce more 
reliable fits that take into account the intrinsic scatter of 
the distribution: 



logE'peak = (4.40 ± 0.509) - (1.31 ± 0.148)r^ 



(3) 



log_Bpoak = (3.05 ± 0.359) - (3.79 ± 0.554) logr^'^^(4) 

Using the values of -Epcak from either a CPL or Band 
function fit is a valid assumption, as shown in Figure|4l This 
figure shows the values of iSpcak derived from both models 
together with the best-fit line, which has a slope very close 
to 1. 

Next we look at all the bursts in the sample, including 
those whose BAT spectra are fit with a CPL model. Since the 
relations above deal with the PL slope of the BAT spectrum, 
we forced the remaining CPL bursts to be fit with a PL 
model to see if the relationships changed in any way (Fig 
[5}. The results are similar, with the slopes becoming slightly 
steeper than the previous values: 

logE'peak = (4.34 ± 0.475) - (1.32 ± 0.129)r^^'^ (5) 

log_Bpoak = (3.00 ± 0.332) - (3.88 ± 0.469) logr^^'^(6) 

Lastly, Figure |6] shows the relationship between p^^^ 
and the joint CPL/Band fit spectral slope a. The data all 
lie below the reference line of a = P^^^. This is understand- 
able, since the effective photon index p^^"^ has to compen- 
sate for the curved spectrum around iJpoak and is, therefore, 
generally steeper. The scatter in the relative steepness is 
related to the location of Ep^nk- The larger the -Epeak, the 
closer ps^"^ to a. If -Epcak is relatively low, the BAT-band 
spectra are severely bent and the effective Pbat could be 
very different from a. 

3.3 Amati Relation 

A sub-sample of 25 GRBs in GBM-BAT joint sample have 
measured redshifts. With this value we can calculate the 
isotropic energy released which, combined with the rest- 
frame peak spectral energy, -Epeak. rest, can be used to test th e 
empirical Amati relation (jAmati et al. 200^ : I Amati 2003h . 
Since our spectral fits are derived from time-integrated spec- 
tra, we can derive the mean energy fiux F during the ex- 
tracted time scale At (which is essentially Tgo) in a desired 
energy band, which we uniformly adopt as the rest-frame 1- 
10* keV energy flux. We worked on only bursts whose time 
integrated spectra are fit by a Band function of CPL, there- 
fore three bursts that are fit with a simple PL in the z-known 



sub-sample are excluded. We then derived the broad-band 
isotropic 7-ray energy -Eiso ~ 4Tvdj^FAt/{l + z), where At 
is the time interval over which the spectral fit was derived, 
d_t is the luminosity distance, and the factor of {1 + z) is 
the correction for cosmological time dilation. If the spec- 
trum is best described by a CPL then the best-fit pre-break 
slope (a) is used with the assumption of a typical /? of -2.3 
to construct a Band function. Otherwise, the best-fit Band 
function parameters are used. In Figure[7]we present the dif- 
ferences between deriving the flux from a CPL as opposed 
to a Band function. The calculated fiuence (and therefore 
Eiso) is systematically underestimated by a mean factor of 
1.4 and a maximum of 2 for this sample. For the bursts that 
are best fit with a CPL we also refit the spectra with a Band 
function and use those spectral parameters, even though this 
model has a larger x^y to derive -Eiso- These results are also 
included in Figure [7] and show that the bursts tend to follow 
a similar, although not identical, distribution as the CPL 
fits. 

We fit the Amati relation sample with two different 
methods: ordinary least squares and ordinary least squares 
bisector fit, the latter an accepted method to fit linear mod- 
els while taking into account some of the scatter of the data 
set. We provide both fits for comparison, although the bi- 
sector method is likely a more accurate depiction of the 
true relation. For the joint sample these methods give dif- 
fering solutions that we explore below. First, using the ordi- 
nary least squares method, our results are generally consis- 



in the literature (LA 


tnati et al. 20081: IZhane et al. 200S : 


Ghirlanda et al. 200^ 
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malization (see Figure |8}. Using the bisector fit we find 
that the slope steepens significantly to 0.76, with a lower 
intercept. This steepening, however, is due to a selection 
effect in the data owing to a lack of bursts above an en- 
ergy of about 10^* erg. The highest energy bursts reported 
by Fermi (e.g. 080916C, 090323) are not simultaneously de- 
tected by Swift and are not included in our sample. Adding 
in the three highest energy Fermi bursts by hand to simu- 
late a more unbiased sample (data from D. Gruber, private 
communication), the relation drops down to 0.57, which is 
fully consistent with previous estimates and serves as a con- 
sistency check for our sample. In general, we derive a larger 
scatter in the Amati relation than the sample presented in 
works such as Amati et al. (2008, 2010). The results are 
consistent, however, with another independent analysis of 
the GBM bursts (Gruber et al. 2011). The mean and me- 
dian of the -Eiso distribution of our sample are 8.4x10^^ 
erg and 6.07 x 10^^ erg, respectively. Figure |S] shows the 
Bpeak.rest — Eiso relationship for the 20 bursts of the joint 
sample with known redshift with the linear bisector fit as 
well as an approximation of the fits by Gruber et al. (2011; 
blue/solid) and Amati et al. (2008; red/dashed) with their 
respective 2a regions. The outlier of the distribution, having 
an Epeak.rost of about 8000 keV, is GRB 090510. It is the only 
redshift-known short burst in the sample. Previous works 
have shown that short {Tgo < 2 sec) bursts follow a different 
relationship and are n ot included in the slope determination 
of the longer bursts jAmati et al. 20081: IZhang et al. 20091 : 
iGhirlanda et al. 20091 : Umati 2010l ). 
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4 COMPARISON OF BURST DURATIONS IN 
THE BAT AND GBM BANDS 

We calculate the Tgo values observed in the GBM band for 
the b ursts in our sam ple by using a Bayesian Block method 
l|Scarg le 1998; Richar dson et al. 1 996). The lightcurves we 
use are in 64 ms bins and are extracted from the brightest 
Nal detector. With the Bayesian Block method we derive 
Ts and Tgs, the epochs when 5% and 95% of the total flu- 
ence were registered, respectively. We then derive Tgo by 
Tqo = Tgs — Ts with the errors of Ts and T95 estimated with 
the bootstrap method and assuming that the observed errors 
in the lightcurves are log-normal. We generate 1000 mimic 
lightcurves from the error distributions and calculate their 
T90 for a given burst. We make a Gaussian fit to the Tgo 
distributions then derive the value of Tgo and its Icr error. 
Tgo il l the BAT band is ta ken from the Second BAT Cat- 
alog (jSakamoto et al. 20111 ) (bursts before 2010) and GCN 
reports (bursts after 2010; see Table 3 for specific references). 

A comparison of the T90 in the BAT and GBM bands 
along with their distributions is shown in Figure It is 
found that the distributions are still bimodal. The Tgo mea- 
sured in the BAT band for long GRBs is usually larger than 
that in the GBM band. However, the T90 measured in the 
BAT band for short GRBs is even shorter than that in the 
GBM band. Since BAT is more sensitive to softer emission 
than GBM, this fact is consistent with the soft /hard na- 
ture of long/short GRBs, respectively. In particular, fainter, 
softer emission can be picked up by BAT but not by GBM, 
so that r9o(BAT) can be much longer than T9o(GBM) in 
some bursts. Conversely, since short GRBs are typically 
hard, BAT may not be sensitive enough to pick up some 
emission episodes that GBM can. As a result, Tgo of short 
GRBs can be longer in the GBM band than in the BAT 
band. 

The difference in T90 for different detectors may also 
affect the derivation of the isotropic energy. This is because 
some long bursts having shorter Tgo in the GBM band than 
BAT band may be sampling only the hardest and bright- 
est part of the spectrum. For the small sample of six bursts 
above the Amati et al. 2g region in Fig. |8] we find that all 
have longer durations in the BAT bands with half of them 
significantly longer (difference > 30%), which may indicate 
the existence of this spectrum sampling effect. This sug- 
gests that in order to fully characterize the durations, and 
hence, the physical origin of the bursts, multi-band data and 
a careful energy-dependent temporal analysis are essential 
(see Qin et al. 2012). 



5 CONCLUSIONS AND DISCUSSION 

Using the data available from Fermi and Swift we performed 
an analysis of the jointly detected GRB sample. The spectral 
analysis indicates that fitting the joint spectra shows con- 
sistency with values reported in the literature base d solely 
on GBM data (|Nava et al. 201ll : lBissaldi et al. 201ll ). and is 
able to reproduce a variety of empirical relations presented 
in the literature, including the V^^^ — -Epcak and Amati re- 
lations. The joint fits do not, in general, give substantially 
different values for the spectral parameters from the GBM- 
only analyses. 



The updated sample of peak energies shows agreement 
with the relation between the PL slope of the BAT spectrum 
and -Epeak derived previously by Sakamoto et al. (2009) and 
Zhang et al. (2007b) and indicates the relation's ability to 
estimate the value of the observed -Epcak from a fit to the 
rather narrow Swift band. The relation shows a fair amount 
of scatter, however, and we caution against using it for ro- 
bust measurements of -Epeak, especially when compared to 
spectra from Fermi or other missions that often contain the 
spectral peak within the energy range of the detectors. It 
may be suited, for example, to give a rough estimate of the 
values of -Epeak for Swift bursts that have no additional ob- 
servations. 

The redshift known subsample of bursts reproduces the 
general relationship in the -Eiso — Speak, rest plane and fur- 
ther confirms the conclusions of Gruber et al. (2011) that 
show generally harder spectra and a larger scatter around 
the relation than previously believed. Li (2007) proposed the 
argument that the -Eiao — Speak, rest relation might evolve with 
redshift. We believe that the sample of bursts with known 
redshift, especially the subsample of joint bursts, is still too 
small to test such a claim. The necessary binning to test the 
hypothesis would thin out the sample sufficiently to make 
firm conclusions tenuous at best. 

Finally, the interesting discrepancy between the Tgo de- 
rived in the GBM and BAT detectors confirms the soft/hard 
nature of long/short GRBs, respectively, and suggests the 
need for a detailed energy-dependent temporal analysis of 
GRBs. 
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Table 1. Joint sample bursts whose time-integrated spectra are best fit by a simple powerlaw model. F is the slope of the 
powerlaw and At the time interval used for the fit. Spectra were fit using data in the energy range of 15-150 keV(BAT)-|-8 

keV-40 MeV (GBM). Reported errors are to the 90% level. 



GRB GBM ID At r X^/dof 

name (s) 



080905B 


080905705 


26.624 


1 62"-''" 


526/530 


080928A 


080928628 


18.432 


^•^'-0.07 


537/529 


090422A 


090422150 


13.312 


9 1 90.24 
^■-^^-0.21 


527/532 


090518A 


090518080 


7.168 


-^•''•^-0.07 


402/421 


090927A 


090927422 


9.216 


1 9302 


689/587 


100619A 


100619015 


102.404 


1 S3OO5 


350/413 


100727A 


100727238 


18.432 


1 qio.og 

-^■''-^-0.09 


545/533 
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Table 2. Bursts whose time-integrated spectra are best fit by either a cutoff powerlaw or Band function model. Spectra 
were fit using data in the 15-150 keV (BAT)-|-8 keV-40 MeV (GBM) energy range. Errors are to the 90% level. 



GRB 


GBM ID 


S])cctral 


At'' 




6 


-^pcak 


xVdof 


uaiiic 




Aludcl" 


(^) 






(k<-\-) 




080714A 


080714745 


GPL 


28.672 


-^■^^-0.15 




155i\i 


401/395 


080725A 


080725435 


CPL 


36.864 


1 1 o0.07 
^•^^-0.08 




3098557 


456/408 


080727C 


080727964 


GPL 


77.825 


J^--'-'J_0.09 




19li'28 


711/526 


080804A 


080804972 


CPL 


26.62 


0.560i2i3 




2001023 


465/408 


08U810A 


080810549 


CPL 


69.633 


1 oqO 06 

-'-■^'^-0.07 




r:eyi269 

5641%o 


581/528 


080905A 


080905499 


GPL 


1.152 


0.25'_!:;fg 




586i\^g2 


283/267 


080916A 


080916406 


GPL 


69.633 


l-24«o.n 




1343^9 


402/410 


081008A 


081008832 


GPL 


60.417 


1-19-0.11 




503!59g7 


591/526 


081012 A 


081012549 


CPL 


18.432 






256LI4S 


506/526 


081024A 


081024245 


CPL 


0.256 


1-05_Q 33 




2170202^074 


104/103 


081102A 


081102739 


GPL 


50.177 




— 


80«„ 


570/529 


081109A 


081109293 


GPL 


40.961 


l-44^-ai5 




133i»3i 


548/530 


081121A 


081121858 


Band 


21.501 


^■^^-0.18 


2-16^o'^l4 


175?.^25 


398/391 


08 11 26 A 


081126899 


GPL 


51.171 






249^347 


425/410 


081221A 


081221681 


Band 


39.422 


ol 


3-14^0.28 


slis 


568/410 


081222A 


081222204 


Band 


14.336 




9 qqO.33 
^.oy_Q 22 


140^13 


460/371 


081226A 


081226044 


CPL 


0.384 








121/120 


090102A 


090102122 


GPL 


37.888 


0.98'_^!!4,, 






499/410 


090113A 


090113778 


GPL 


12.288 


-4 ^ in 

1 oqu.iy 
^•^^-0.23 




138L°3°9 


346/348 


090129A 


090129880 


GPL 


15.36 






163?.^ 


343/376 


090423A 


090423330 


GPL 


14.336 


0.870-3''^3 


_ 


52^6 


517/506 


090424A 


090424592 


Band 


19.71 


0.95002^2 


9 SqO.26 
^•oy_0.17 


155^5 


844/525 


090509A 


090509215 


GPL 


26.624 


0-9°,l23 




221 "Is 


572/527 


090510A 


090510016 


GPL 


0.896 


0-84'^ao5 




4482f5rfsi 


329/304 


090519A 


090519881 


GPL 


18.432 


f) Ki 0.28 




296 ^"^^3 


532/518 


090531B 


090531775 


CPL 


3.072 


0.930i5i6 




191lL"|'/3 


371/347 


090621B 


090621922 


CPL 


0.128 


0-47°0.34 




778 26 
^' '-230 


57.2/77 


090708A 


090708152 


CPL 


14.336 






54L\ 


460/496 


090709B 


090709630 


CPL 


20.48 






1211^20 


507/530 


090712A 


090712160 


GPL 


53.248 


l-06-ai2 




Q Q 

5«»_204 


492/529 


090813A 


090813174 


GPL 


9.725 


l-51-"i4 




lOO'l^s 


344/338 


090904B 


090904058 


GPL 


74.663 


1 oq().(J8 
-'-^*-0.08 


— 


139L'\4 


554/529 


090912A 


090912660 


CPL 


123.902 


sO.22 
^•°-0.24 




72?. 7 


536/530 


090926B 


090926914 


CPL 


64.513 


90.14 

"■^-0.15 


— 


83^4 


520/529 


091020A 


091020900 


CPL 


28.672 


1 Qfi0.08 

1.0D_Q Qg 




2m%^ 


382/412 


091026A 


091026550 


GPL 


15.36 


1 Q10.15 


— 


n'yn235 
^'^-99 


410/378 


091102A 


091102607 


GPL 


9.216 


l-03°ai7 


— 


439164^7 


514/521 


091127A 


091127976 


Band 


11.006 


1-21°J:'i5 


9 9c0.04 
^•^O_0.04 


35^3 


466/366 


091208B 


091208410 


CPL 


13.312 


1.0D_Q Q9 




118-15 


486/375 


091221A 


091221870 


CPL 


35.841 


■^■^'-0.07 


— 


275%^ 


721/529 


lOOlllA 


100111176 


CPL 


16.374 






203i_2222 


706/649 


100117A 


100117879 


GPL 


0.192 


n 1 4.0-39 


— 


380i'\4^2 


93.1/86 


100206A 


100206563 


GPL 


0.256 


n fi4.o.2 

^•'^^-0.21 




763?.f69 


109/126 


100504A 


100504806 


GPL 


28.672 


^•^^^-0.24 




85-11 


582/528 


100522A 


100522157 


GPL 


41.985 


1 c^O 22 
1-54^0 18 




59 9 


393/412 


100615A 


100615083 


Band 


40.574 


910-29 
"■^ -0.38 


2 070- 20 

-0.09 


50-\o 


437/410 


100625A 


100625773 


CPL 


0.32 


69011 




^^'^-131 


160/143 


100704A 


100704149 


CPL 


23.552 


910 09 




1873323 


584/528 


100728A 


100728095 


GPL 


191.487 


n «o.03 




313L«i7 


538/411 


100728B 


100728439 


GPL 


10.237 


0.930182 




10522j,^ 


537/501 


100802A 


100802240 


GPL 


41.985 


0-57"_-§^4i 




832^^2 


684/650 


100814A 


100814160 


GPL 


158.722 


1.060_i\2 




14733^^ 


566/526 


100816A 


100816026 


CPL 


4.542 


50.11 




13913^1 


437/364 


100906A 


100906576 


Band 


121.347 


1 440.11 


1 QqO.23 


10635,4 


589/526 


100924A 


100924165 


GPL 


8.061 


^•^"-0.17 




1 72*'*^ 


472/437 


101008A 


101008697 


GPL 


14.336 


1 24"'l'' 
-^•^^-0.22 




70l5^,r 


422/410 


lOlOllA 


101011707 


GPL 


27.645 


1 0402 




389000^^ 


520/524 
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Table 2 (cont'd) 



GRB 


GBM ID 


Spectral 


At'' 


a 


/3 -Epoak 




name 




Model" 


(s) 




(keV) 




101024A 


101024486 


CPL 


13.312 




48^5 


543/504 


101213A 


101213451 


GPL 


47.105 


1.260^1.12 


423i\043 


380/411 


101219B 


101219686 


CPL 


51.201 


0.08^^39^^ 


72^3 


738/650 


110102A 


110102788 


CPL 


148.486 




47969^1,^ 


403/411 


110106B 


110106893 


CPL 


35.84 


1 4().18 




509/528 


110119A 


110119931 


CPL 


68.609 


1 Infill 


224Ll°s 


555/527 


110201A 


110201399 


CPL 


12.288 




7306«j 


612/595 


110213A 


110213220 


CPL 


43.009 


1.560-09i 




411/410 


110318A 


110318552 


Band 


20.48 


^■82%\^ 


n qO.4 Qo13 
'^•'^-0.18 ^^-11 


446/394 


110402 A 


110402009 


CPL 


38.526 




9143_248g 


618/649 


110412A 


110412315 


GPL 


20.48 


0.640J4g 


84^2, 


532/528 



''GPL = CutoflF powerlaw or powerlaw + exponential. 100 keV normalization energy. 

''Duration of time bin for the evaluation of the time-integrated spectra. This is determined 
by the duration of the signal of the GBM bursts in RMfit 

■^Several bursts in the table (080727G, 081221, 081222, 090424, 091127, 091208B, 100728A) 
have values larger than about 1.3. These bursts are usually long (20 s > At > 200s) and/or 
have multi-peak structure, and show significant spectral evolution over the duration of the burst. 
Doing a simple time-dependent analysis show the fits improve but the spectral evolution likely 
drives the larger residuals for the full time interval. 
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Table 3. GBM and BAT Tgo values in the 8-1000 keV and 15-150 keV energy ranges, respectively. GEM values are derived 
with the method described in §2.2 and BAT figures are referenced from the literature (see reference column for appropriate 

citations). 



GRB 

name 


GBM ID 


GBM Tgo 

(s) 


BAT Tgo 

(s) 


Reference* 


080714 


GRB080714745 


6.34±0.36 


25.81 


bat2 


080725 


GRB080725435 


22.21±0.23 


92.74 


bat2 


080727C 


GRB080727964 


35.55±0.59 


77.61 


bat2 


080804 


GRB080804972 


73.41±0.93 


37.19 


bat2 


080810 


GRB080810549 


49.34±0.63 


107.67 


bat2 


080905A 


GRB080905499 


1.06±0.30 


1.02 


bat2 


080905B 


GRB080905705 


192.8±1.15 


101.62 


bat2 


080916A 


GRB080916406 


44.26±0.72 


61.35 


bat2 


080928 


GRB080928628 


24.54±0.44 


233.66 


bat2 


081008 


GRB081008832 


175.2±1.16 


179.52 


bat2 


081012 


GRB081012549 


12.8±0.56 


25.2 


bat2 


081024A 


GRB081024245 


0.13±0.18 


1.82 


bat2 


081025 


GRB081025349 


23.62±0.49 


22.78 


bat2 


081101 


GRB081101491 


0.54±0.39 


0.18 


bat2 


081102 


GRB081102739 


29.47±0.59 


48.19 


bat2 


081109A 


GRB081109293 


27.46±0.65 


221.49 


bat2 


081121 


GRB081121858 


17.98±0.42 


17.67 


bat2 


081126 


GRB081126899 


35.36±0.46 


57.65 


bat2 


081221 


GRB081221681 


45.82±0.99 


33.91 


bat2 


081222 


GRB081222204 


26.75±0.53 


33 


bat2 


081226A 


GRB081226044 


0.48±0.20 


0.44 


bat2 


090102 


GRB090102122 


29.02±0.54 


29.3 


bat2 


090113 


GRB090113778 


8.58±0.29 


9.1 


bat2 


090129 


GRB090129880 


14.02±0.23 


17.66 


bat2 


090422 


GRB090422150 


- 


8.47 


bat2 


090423 


GRB090423330 


12.35±0.50 


9.77 


bat2 


090424 


GRB090124592 


45.79±0.61 


49.47 


bat2 


090509 


GRB090509215 


261.18±1.08 


336.38 


bat2 


090510 


GRB090510016 


0.38±0.045 


5.66 


bat2 


090516A 


GRB090516353 


85.28±0.54 


208 


bat2 


090518 


GRB090518080 


7.97±0.44 


85.82 


bat2 


090519 


GRB090519881 


42.88±0.57 


58.26 


bat2 


090531B 


GRB090531775 


2.016±0.39 


55 


bat2 


090618 


GRB090618353 


130.24±1.05 


113.34 


bat2 


090621B 


GRB090621922 


0.29±0.16 


0.14 


bat2 


090708 


GRB090708152 


12.48±0.39 


8.7 


bat2 


090709B 


GRB090709630 


11.14±0.32 


27.02 


bat2 


090712 


GRB090712160 


31.68±0.62 


186.68 


bat2 


090813 


GRB090813174 


8.96±0.55 


7.14 


bat2 


090904B 


GRB090904058 


52.35±0.59 


64 


bat2 


090912 


GRB090912660 


126.37±0.54 


135.52 


bat2 


090926B 


GRB090926914 


41.34±0.62 


99.28 


bat2 


090927 


GRB090927422 


3.2±0.32 


2.16 


bat2 


091020 


GRB091020900 


44.67±0.65 


38.92 


bat2 


091024 


GRB091024372 


48.26±0.62 


112.28 


bat2 


091026 


GRB091026550 


15.74±0.46 


174 


bat2 


091102 


GRB091102607 


7.58±0.46 


6.65 


bat2 


091112 


GRB091112737 


50.85±0.82 


19.59 


bat2 


091127 


GRB091127976 


9.15±0.26 


7.42 


bat2 


091208B 


GRB091208410 


11.39±0.14 


14.8 


bat2 


091221 


GRB091221870 


30.432±0.62 


68.49 


bat2 


lOOlllA 


GRB100111176 


11.01±0.46 


12.9 


10322 


100117A 


GRB100117879 


0.51±0.19 


0.3 


10338 


100206A 


GRB 100206563 


0.19±0.13 


0.12 


10379 


100212 


GRB100212588 


8.10±0.46 


136 


10404 


100413A 


GRB100413732 


78.11±1.07 


191 


10600 


100427A 


GRB100427356 


10.43±0.41 
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Table 3 (cont'd) 



GRB 


GBM ID 


GBM Tgo 


BAT Tgo 


Reference" 


name 




(s) 


(s) 




100504A 


GRB100504806 


24.51±0.73 


97.3 


10716 


100522A 


GRB100522157 


37.38±0.59 


35.3 


10788 


100615A 


GRB100615083 


36.42±0.59 


39 


10850 


100619A 


GRB100619015 


92.61±0.34 


97.5 


10864 


100625A 


GRB100625773 


0.38±0.14 


0.33 


10891 


100704A 


GRB100704149 


11.81±0.54 


197.5 


10932 


100727A 


GRB 100727238 


25.92±0.47 


84 


11001 


100728A 


GRB 100728095 


159.97±0.76 


198.5 


11018 


100728B 


GRB100728439 


9.34±0.46 


12.1 


11023 


100802A 


GRB100802240 


132.26±0.84 


487 


11035 


100814A 


GRB100814160 


25.41±0.29 


174.5 


11094 


100816A 


GRB100816026 


2.24±0.23 


2.9 


11111 


100906A 


GRB100906576 


115.56±0.68 


114.4 


11234 


100924A 


GRB100924165 


11.87±0.43 


96 


11296 


101008A 


GRB101008697 


7.2±0.30 


104 


11327 


lOlOllA 


GRB101011707 


38.50±0.62 


71.5 


11332 


101023A 


GRB101023951 


66.94±0.55 


80.8 


11367 


101024A 


GRB101024486 


20.51±0.54 


18.7 


11374 


101129A 


GRB101129652 


0.544±0.17 


- 


- 


101201A 


GRB101201418 


83.62±0.68 


- 


- 


101213A 


GRB101213451 


38.24±0.56 


135 


11453 


101219B 


GRB101219686 


54.21±0.75 


34 


11475 


101224A 


GRB101224227 


0.56±0.32 


0.2 


11486 


110102A 


GRB110102788 


133.76±0.34 


264 


11511 


110106B 


GRB110106893 


22.78±0.62 


24.8 


11533 


110119A 


GRB110119931 


59.87±0.59 


208 


11584 


110128A 


GRB110128073 




30.7 


11614 


110201A 


GRB110201399 


11.712±0.51 


13 


11624 


110207A 


GRB110207470 


39.01±0.52 


80.3 


11664 


110213A 


GRB110213220 


33.12±0.46 


48 


11714 


110318A 


GRB110318552 


14.66±0.41 


16 


11802 


110319B 


GRB110319815 


13.95±0.52 


14.5 


11818 


110402A 


GRB110402009 


35.19±0.53 


60.9 


11866 


110412A 


GRB110412315 


17.4±0.53 


23.4 


11929 


110420B 


GRB110420946 


0.74±0.25 


0.084 


11946 



'"References for BAT Tgo: bat2 - Second BAT Catalog 
l lSakamoto et al. 201lh . Other numbers in this column correspond to 
GCN reports. See bibliography for full record information. 
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Figure 1. Distributions of ii^neak a for the best fit joint sample bursts. Both are generally consistent with the literature (e.g. 
IPreece et al. 200ol : IZhang et al. 201ll : iNava et al. 201lh . The solid (black) and dashed (red) histograms are for long and short GRBs, 
respectively. 




E , Joint [keV] 



Figure 2. The relationship between the best fit Speak for the joint sample and the GBM spectra. The line superimposed is the line with 
slope 1. This shows that adding the BAT data into the fit does not significantly change the values of i?pcak- 
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Figure 3. The V^^^ - 

^pcak relation in both log-log (a) and log-linear (b) space created from 53 best-fit spectra. 15 bursts were removed 
from the original sample of 75 because the BAT best-fit model was a CPL and not a simple PL. An additional seven bursts were removed 
because the joint spectrum is described by a simple PL. As shown below, the difference in Bpeak between the CPL and Band model fit for 
the joint spectra is negligible, so the best-fit model value was used. The lines superimposed on the plots are the best-fit lines generated 
using an ordinary least squares fit (red, solid) and a least squared linear bisector fit (blue, dash) that considers some of the scatter in 
the relation. The latter procedure's parameters are summarized in the text. 
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Figure 4. The relationship between the values joint Epoak for a- CPL and Band model for all 75 GRBs in our sample. The best fit line 
(solid, black) has a slope of nearly 1, as shown by the reference line (dashed, red), with very few outliers from this relation. Such small 
changes in the value are negligible in the global analysis of spectral properties and we, therefore, always take the best-fit model value of 
Speak- Some Band function fits are not well constrained and are reflected in the large horizontal error bars for some bursts. 
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Figure 5. The r^^''^ — -Epcak relation in both log-log (a) and log-linear (b) space created from all the bursts in our sample. We refit 
the Swift bursts that are described as a CPL with a simple PL model and use this value of F, regardless of the increase in the x^- The 
lines superimposed on the plots are the best-fit lines generated using an ordinary least squares fit (red, solid) and a least squared linear 
bisector fit (blue, dash) that considers some of the scatter in the relation. The latter procedure's parameters are summarized in the text. 
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Figure 6. Relationship between the BAT PL spectral slope and the joint fit CPL/Band a for short (red, triangle) and long (black, 

square) bursts. The line superimposed is a reference of slope = 1. 
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Figure 7. Alternative methods for calculating the relationship between Eiso jcst-Epcak.rcst- In both panels, the black (square) icons and 
the solid black line are identical to those presented in Figure 8. The dashed line represents the best linear bisector fit to the second data 
set for reference. Short burst 090510 is labelled but not included in the calculations, (a) Rest frame fluxes are calculated assuming a 
CPL spectrum (red/circle) instead of a Band function (black/square, identical to Fig 8). This systematically underestimates the isotropic 
energy by a mean factor of about 1.4. (b) Rest frame fluxes calculated using Band function parameters, even if a CPL gives a better 
fit (i.e. lower x'^; blue/triangle). If /3 cannot be constrained a typical value of -2.3 is assumed. This gives a similar distribution to that 
produced in Figure 8 (black/square). 
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Figure 8. Relationship between the isotropic equivalent energy release and the peak energy of the time-integrated spectrum of the joint 
sample bursts. Isotropic energies are derived using a Band function spectrum. If the spectrum is best described by a CPL then the best-fit 
pre-break slope (a) is used with the assumption of a typical /3 of -2.3 to construct a Band function. Superimposed are the best ordinary 
least squares bisector fit of the joint sample (black/dash-dot, slope=0.76) and the approximate best fit lines and 2a regions reported in 
Amati et al. (2008; red/dashed; slope=0.57) and Gruber et al. (2011; blue/solid; slope=0.52). The shift to steeper slopes is produced 
from a lack of bursts above Eiso ~ 10^** erg, which is attributed to a selection eff'ect of the sample. These bursts are the highest energy 
bursts detected and are Fermi LAT triggers and therefore do not show up in our sample of joint GBM-BAT bursts. If these highest 
energy bursts are introduced to create a more complete sample, we show consistency with both slope values presented in the literature. 
We also confirm the larger scatter and slightly harder spectra presented by Gruber et al. (2011). Differences in the absolute values of 
^poak s-re likely due to the derivation of the parameter using the joint BAT-GBM data. The obvious outlier in the distribution is GRB 
090510 (labelled), the only short burst (Tgo < 2 sec) in the redshift-known subsample. This burst is not used in calculating the best-fit 
line. Previous works (Amati et al. 2008; Zhang et al. 2009; Ghirlanda et al. 2009; Amati et al. 2008) have shown that these bursts follow 
a different relationship in the Epeak.rest ~ ^-iso plane and are not considered in the fitting of the slope of the relation. 
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Figure 9. Tgo measures for the joint sample bursts as derived witii the BAT (15-150 keV) and GBM (8-1000 kcV)data and corresponding 
distributions. The diagonal solid line is the line of slope = 1, while the dashed lines indicate a duration of 2 sec. Note the typically longer 
durations in the Swift band for the long and generally softer bursts while the trend reverses for the shorter, typically softer, bursts. 



